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Abstract 

Starting from the basic prescriptions of the Tsallis' nonextensive thermostatistics, 
i.e. generalized entropy and normalized (/-expectation values, we study the relativis- 
tic nonextensive thermodynamics and derive a Boltzmann transport equation that 
implies the validity of the i?-theorem where a local nonextensive four-entropy den- 
sity is considered. Macroscopic thermodynamic functions and the equation of state 
for a perfect g derived at the equilibrium. 

PACS: 05.20.Dd, 05.70.Ln, 05.90.+m, 25.75.-q 



1 Introduction 



Recently, there is an increasing evidence that the generalized nonextensive 
statistical mechanics, proposed by Tsallis [1], can be considered as the more 
appropriate basis of a theoretical framework to deal with physical phenom- 
ena where long-range interactions, long-range microscopic memories and/or 
fractal space-time constraints are present (cf. [1] for details). A considerable 
variety of physical applications involve a quantitative agreement between ex- 
perimental data and theoretical models based on Tsallis' thermostatistics [2]. 
In particular there is a growing interest to high energy physics applications 
of nonextensive statistics. Several authors outline the possibility that experi- 
mental observations in relativistic heavy-ion collisions can reflect nonextensive 
features during the early stage of the collisions and the thermalization evolu- 
tion of the system [3-7] . 

The basic aim of this letter is to study the nonextensive statistical mechanics 
formalism in the relativistic regime and to investigate, through an appropri- 
ate relativistic Boltzmann equation, the non-equilibrium and the equilibrium 
thermodynamics relations. 
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2 Basic assumptions in nonextensive thermostatistics 



Let us briefly review some basic assumptions of the nonextensive thermostatis- 
tics that will be useful in view of the relativistic extension. 

Starting point of the Tsallis' generalization of the Boltzmann-Gibbs statistical 
mechanics is the introduction of a g-deformed entropy functional defined, in 
a phase space system, as [1] 

S q = -k B f dttp q ln q p, (1) 



where k B is the Boltzmann constant, p = p(x, v) is the phase space probability 
distribution, dfl stands for the corresponding phase space volume element and 
ln q x = (rr 1 " 9 — 1) / (1— q) is, for x > 0, the g-deformed logarithmic function. For 
the real parameter q — > 1, Eq.(l) reduces to the standard Boltzmann-Gibbs 
entropy functional. 

In the equilibrium canonical ensemble, under the constraints imposed by the 
probability normalization 



J dttp = 1 , (2) 



and the normalized g-mean expectation value of the energy [8] 

SdnjfiH(x,v) 

the maximum entropy principle gives the probability distribution [8,9] 

P{x,v) = — — , (4) 



where 

f(x, v) = [1 - (1 - q)P(H(x, v) - (H),)] 1 ^ (5) 



and 



Z q = J dQf(x,v) . (6) 



Let us note that, depending from the extremization procedure, in Ref. [8] the 
above factor (5 is only proportional to the Lagrange multiplier, because of the 
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probability distribution is self-referential, while, in Ref.[9], (5 is actually the 
Lagrange multiplier associated to the energy constraint. 

By using the definitions in Eqs. (4) and (6) into the relation Z l q ~ q = J dVt p q 
[9], the following identity holds [10] 

J<mf(x,v) = Jdnr(x,v), (7) 

and the normalized g-mean expectation value for a physical observable A(x, v) 
can be expressed as 

Therefore, the probability distribution and the g-mean value of an observable 
have the same normalization factor Z q , as in the extensive statistical mechan- 
ics. Such a non-trivial property does not depend on the equilibrium frame 
but it comes from the normalization condition and holds at any time (if we 
require that the transport equation conserves the probability normalization or 
the number of particles). This observation will play a crucial role in the correct 
formulation of the relativistic Boltzmann equation and in the definition of the 
thermodynamic variables. 



3 Relativistic kinetic theory 

On the basis of the above prescriptions, we are able to progress in the formula- 
tion of the relativistic nonextensive statistical mechanics. Let us start defining 
the basic macroscopic variables in the language of relativistic kinetic theory. 
Because we are going to describe a non-uniform system in the phase space, we 
introduce the particle four-flow as 

N»(x) = ± f^p»f(x,p), (9) 

Zi q J p 

and the energy-momentum flow as 

T^(x) = ^ f^-p^r(x,p), (10) 

Zi q J p 

where we have set h — c — 1, x = x^ = (t, x), p = p^ = (p°,p) and 
p° = y/p 2 + m 2 is the relativistic energy. The four- vector JV = (n,j) con- 
tains the probability density n = n(x) (which is normalized to unity) and 
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the probability flow j = $(x). The energy- momentum tensor contains the nor- 
malized g-mean expectation value of the energy density, as well as the energy 
flow, the momentum and the momentum flow per particle. Its expression fol- 
lows directly from the definition (8); for this reason it is given in terms of 
f q (x,p). 

In order to derive a relativistic Boltzmann equation for a dilute system in 
nonextensive statistical mechanics, we consider the finite volume elements 
A 3 x and A 3 p in the phase space. These volume elements are large enough 
to contain a very large number of particles but also small enough compared 
to the macroscopic dimension of the system. If, in the volume A 3 xA 3 p, is 
contained a representative sample of system, we can assume that Eq.(7) still 
holds in such phase space portion. Then, in the Lorentz frame, the particle 
fraction AN(x,p) in the volume A 3 a;A 3 p can be written as 

AN(x,p) = ^- J J dQp(x,p), (11) 

where N is the total number of particles of the system. With respect to the 
observer's frame of reference, the above expression becomes 

AN(x,p) = ^J J d 3 a^p»nx,p), (12) 

q A 3 o- A a p 

where d 3 a VL is a time-like three-surface element of a plane space-like surface a 
[11]. Requiring that the net flow through the surface A 3 a of the element A 4 x 
vanishes in absence of collisions, we have: p^d fl f g (x,p) = 0. While considering 
collisions between particles, the Boltzmann equation becomes 

p^r(x,p) = C q (x,p), (13) 



where C q (x,p) is the g-deformed collision term that, under the hypothesis that 
only binary collisions occur in the gas, can be expressed as 

n i \ 1 fd 3 pid 3 p' d 3 p\ ( , , , , 
zj Pi P Pi 

-h g [fJi]W(p,p 1 \p',p' 1 )} . (14) 

In the above equation, we have set W(p, Pi\p',p'i) as the transition rate be- 
tween two particle state with initial four-momentum p and p' and a final state 
with four- momenta p\ and p\\ h q [f, fi] is the correlation function related to 
two particles in the same space-time position but with different four-momenta 
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p and pi, respectively. The factorization of h q in two single probability distri- 
butions (uncorrelated particles at the same spatial point) is the celebrated hy- 
pothesis of molecular chaos (Boltzmann's Stosszahlansatz). Thus, the function 
h q defines implicitly a generalized nonextensive molecular chaos hypothesis. 

By assuming the conservation of the energy-momentum in the collisions (i.e. 
p^+p'^ = Pi+p'i) and requiring that the correlation function h q is symmetric 
and always positive (h q [f, fi] = h q [f±, /], h q [f, f\] > 0), it is easy to show that 
collision term satisfies the following property 

m = /^^,p)C ? (x,p) = 0, (15) 

if 

i>(x,p) = a(x) + b II (x)p^ , (16) 
where a(x) and b(x) are arbitrary functions. 

By choosing ip = const, and using the Boltzmann equation (13), then Eq.(15) 
implies 

-Jdnr(x, P ) = o, (17) 

and this is, on account of Eq.(7), nothing else that the conservation of the 
probability normalization Z q . Otherwise, by setting ip = b^, we have from 
Eq.(15) 

d v T^{x) = , (18) 

which implies the energy and the momentum conservation. 

Let us remark that to have conservation of the probability normalization, 
energy and momentum, it is crucial that not only the collision term C q be 
explicitly deformed by means of the function h q , but also the streaming term 
P^dpf 9 - This matter of fact is a directly consequence of the nonextensive sta- 
tistical prescription of the normalized g-mean expectation value (8) and is not 
taken into account in the non-relativistic formulation of Ref.[12]. 
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4 Local //-theorem 



The relativistic local if-theorem states that the entropy production cr q (x) = 
dpSgfe) at any space-time point is never negative. 

Assuming the validity of the Tsallis entropy (1), it appears natural to introduce 
the nonextensive four- flow entropy Sff(x) as follows 

W = ~K I ^-jff q (x,p)[\n q f(x,p) - 1] . (19) 
J p 

On the basis of the above equation the entropy production can be written as 

a q (x) = -k B f^ ln 9 / jfdpf = -k B F[]n q /] , (20) 
J p 

where the second identity follows from the Boltzmann equation (13) and the 
definition of F{ip] in Eq.(15). After simple manipulations, Eq.(20) can be 
rewritten as 

k B r d 3 pi d 3 pi d 3 p' d 3 pi' 
<r q W = YJ p o p o 75~^o- 

{h q [f'J'i]W(p',p' 1 \p,p 1 )-h q [fJ 1 }W(p,p 1 \p',p' 1 )}. (21) 

By assuming the detailed-balance property W(p, Pi\p',p'i) = W(p' ,p\\p,Pi), 
we have that the entropy production is always an increasing function, if q > 
and if the function h q satisfies the general condition 

h q [f,f 1 ] = h q [\n q f + \n q f 1 ], (22) 
in addition to be symmetric and always positive. 

Because the nonextensive formalism reduces to the standard Boltzmann ki- 
netic formulation for q — > 1, it appears natural to postulate the g-generalized 
Boltzmann molecular chaos hypothesis as 

fc s [/,/i]=e ff (ln,/ + ln,/ 1 ), (23) 

where we have introduced the Tsallis g-exponential function 

e q (x) = [1 + (1 - q)x} 1 ^ , (24) 



(ln g f + ln 9 / / 1 -ln 9 / + ln g / 1 ) x 
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which satisfies the properties: e q (\n q x) = x and e q (x) ■ e q (y) = e q [x + y + 
(1 — q)xy\. Let us note that a similar expression for the function h q was previ- 
ously introduced in Ref.s[12,13] and a rigorous justification of the validity of 
the ansatz (23) can be found only by means of a microscopic analysis of the 
dynamics of correlations in nonextensive statistics. 



5 Equilibrium and equation of state 



The condition that entropy production vanishes everywhere, together with re- 
quirement that the equilibrium probability distribution f eq must be a solution 



of the transport equation (13) 
ing into account of Eqs.(15), 
occur when ln q f eq = a + b^p 11 
after simple redefinition of the coefficients a and b, the equilibrium probability 
distribution can be written as a Tsallis-like distribution 



, uniquely defines the state of equilibrium. Tak- 
(16) and (20), the condition o q = can only 
By imposing that f eq must satisfy Eq.(13) and 



feq(p) = Y 



1 -(!-?) 



1 1/(1-9) 



k B T 



(25) 



where is the hydrodynamic four- velocity [11] and f eq depends only on the 
momentum in absence of an external field. At this stage, k B T is a free parame- 
ter and only in the derivation of the equation of state it will be identified with 
the physical temperature. Moreover, it is easy to show that f eq is a solution 
of the transport equation (13). 

We are able now to evaluate explicitly all other thermodynamic variables and 
provide a complete macroscopic description of a relativistic system at the 
equilibrium. Let us first calculate the probability density defined as 

n = N^U, = A- / ^ p»U, f eq (p) . (26) 
Zj q j p 



Since n is a scalar, it can be evaluated in the rest frame where £/ M = (1, 0, 0, 0). 
Setting r = p°/k B T and z = m/k B T, the above integral can be written as 



A 

q z 
A °° 



(27) 
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where the last identity has been obtained by a partial integration. Let us 
introduce the g-modified Bessel function of the second kind as follows 

A '»<«- z) = S f /* (t2 - z2) '" ir ' ( e « r )' • (28) 



then, the particle density can be cast into the compact form 

4"7r 

n = — m 2 k B TK 2 (q,z) . (29) 



Similarly we can obtain the other macroscopic thermodynamic variables. Con- 
sidering the decomposition of the energy- momentum tensor [11]: T^ v = e U^U U - 
pA^, where e is the energy density, p the pressure and = g^ v — U^U", 
the equilibrium pressure can be calculated as 

P = -\^ V A, V = J ^p»\ v f e \{p) , (30) 



and can be expressed as 



Am- 

p=^m 2 (k B T) 2 K 2 (q,z). (31) 



Comparing Eq.(29) with Eq.(31), we obtain 

p = nk B T, (32) 

which is the equation of state of a perfect gas if we identify T as the physical 
temperature of the system. An explicit g-dependence does not appear. The 
same equation has been derived in non-relativistic regime and in nonextensive 
scenario [10,14]. 



We proceed now to calculate the energy density e as 

1 f d 3 p i 

7j q 



= T^U,U V = ±1 ^(p»U»Yf? q (p) . (33) 

Zj q J p 



Inserting the previously defined variables r and z and using the definition in 
Eq.(28), we obtain 



4tt 4 
e = — m 



(34) 
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Thus the energy per particle e = e/n is 



e = 3 k B T + m 



K 2 (q,z) ' 



(35) 



which has the same structure of the relativistic expression obtained in the 
framework of the equilibrium Boltzmann-Gibbs statistics [11]. 

In the non-relativistic limit (p< 1) the energy per particle reduces to the 
well-known expression 



Also in this case, no explicit g-dependence is detected. 



6 Conclusion 

The physical motivation of this investigation lies in the strong relevance that 
nonextensive statistics could have in high energy physics. In this letter we have 
studied the thermostatistics of a relativistic system in nonextensive statistics; 
the obtained results can be easily extended to the case where an external force 
is present. 

The prescription of the normalized g-mean expectation values implies a con- 
sistent nonextensive generalization of the macroscopic variables and T^ u . 
On this basis, we have derived a generalized Boltzmann equation where both 
the streaming and the collision terms depend on the deformation parameter q. 
Such a transport equation conserves the probability normalization (or number 
of particles) and is consistent with the energy-momentum conservation laws. 
The collision term contains a generalized expression of the molecular chaos and 
for q > implies the validity of a generalized if-theorem, if the nonextensive 
local four-density entropy (19) is assumed. At the equilibrium, the solution of 
the Boltzmann equation is a relativistic Tsallis-like distribution and the equa- 
tion of state of a classical relativistic gas in nonextensive statistical mechanics 
has the same form as in ordinary Boltzmann-Gibbs frame. 

Finally, nonextensive statistical effects can be detected in connection to micro- 
scopic observables such as particle distribution, correlation functions, fluctu- 
ations of thermodynamical variables but not directly in connection to macro- 
scopic variables, such as temperature or pressure, because the equation of state 
and the macroscopic thermodynamical relations do not explicitly depend on 



e ~ m H — k„T . 
2 B 



(36) 
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the deformation parameter q. 
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